a^ 

CD 

m 

(N 



> 
00 
00 

o 

o\ 
o 

(D 



X 



CERN 97-256 
DESY 97-173 

Non-perturbative 0(a) improvement of the vector current * 

M. GuagncUi'' and R. Sommer'' 

'^CERN, Theory Division, 1211 Geneva 23, Switzerland 
'^DESY-IfH, Platanenallee 6, D-15738 Zeuthen, Germany 

We discuss non-perturbative improvement of the vector current, using the Schrodinger Functional formalism. 
By considering a suitable Ward identity, we compute the improvement coefficient which gives the 0(a) mixing of 
the tensor current with the vector current. 



1. INTRODUCTION 

It is well known that in Wilson's lattice QCD 
discretization errors receive contributions that 
are linear in the lattice spacing. The most obvi- 
ous recipe to reduce them — namely to go to the 
continuum limit by performing numerical simula- 
tions at smaller lattice spacings — represents still 
a hard task, even in the quenched approximation. 

Following Symanzik [|l| , lattice artifacts can be 
removed order by order in a by adding appro- 
priate higher-dimensional operators to the action 
and to the fields whose correlation functions are 
of interest. If one restricts oneself by requiring 
improvement only for on-shell quantities |^ , such 
as particle masses and matrix elements between 
physical states, the structure of the improved ac- 
tion for QCD and of the improved currents is 
rather simple. The general theory was developed 
in detail in ref. In the following we will use 
the same notation without further reference. 

In the quenched approximation, the improved 
action was determined non-perturbatively for 
couplings < go < I 0. 

The renormalized and improved axial-vector 
and vector currents are parametrized as 



(^r)m 
with 



ZA(H-5Aamq)(Ai);i, 
Zv(l + 6vamq)(Vi);, 



ya 



(^l)^ 
(Vi); 

In the above expressions, the coefhcients cx 
and bx are functions of the bare coupling, and 

*Talk given by M. Guagnelli at the International Sym- 
posium on Lattice Field Theory, 21-27 Juli 1997, Edin- 
burgh, Scotland 



are known to 1-loop order of perturbation thc- 
ory[^|j|]. Moreover, the normalizations {Za and 
Zy) and the improvement coefficients ca and by 
have been determined non-perturbatively in the 
quenched approximation for go < 1 As a 

next step, we determine cy, thus completing the 
knowledge of the improved vector current. The 
computation of other improvement coefficients is 
discussed in ||^. 

Knowledge of cy is, for example, required for 
the computation of vector meson decay constants. 
The relative contribution of ad^T^^u to these de- 
cay constants can be, at = 1, as large as 
0.3 X Cy Although the perturbative esti- 

mate 1^, 

cy = -0.01225(1) X glCY + 0{gt), Cp = 4/3,(1) 

suggests an effect of less than 1% on the decay 
constants, our preliminary non-perturbative re- 
sults determine cy to be much larger in magni- 
tude for On — 1- 



2. THE STRATEGY 

Chiral Ward identities relate correlation func- 
tions of axial-vector and vector currents. In the 
0(a) improved theory and for zero quark mass 
these identities can be written in a form which is 
valid up to error terms that are quadratic in a. 
Since the 0(a)-improved axial-vector current as 
well as Zy and Zj^ are known, one can use a par- 
ticular Ward identity in which the only unknown 
is Cy. 

Our starting point is the Ward identity (in the 
continuum) : 



Aa,{x)t'^''%Al{x)Al{y)Q'^) 



(2) 



dR 



2 



2m / d''xe'''>%P''{x)Al{y)Q^) =2i{V^{y)Q^) , 

JR 

where the space-time region R with boundary dR 
contains the point y and is a source located 
outside R. 

We then impose Schrodinger Functional bound- 
ary conditions , choose v = k and specify the 
source, 



(3) 



where C are the quark fields at the boundary 
2^0 = . Finally the region R is specified to be 
the space-time volume between the hyperplanes 
= ti and = ^2- 

On a lattice with finite spacing a, we then ex- 
pect that 

a3 5]e'^''^([(ARX(t2,x)-(ARX(ti,x)]x 

{A^)l{y)Ql) = 

2« ((^R)fc(y)QD + 0(a2), <i < 2/0 < t2 , (4) 

is valid with correction terms of order a? . 

Note that passing from eq. (|^) to eq. (Q) we 
have specialized to to = 0. The reason for this 
is as follows. For finite quark mass, the second 
term on the left hand side of eq. (H) contains zero 
separations x — y. Since it is not an integral over 
an on-shell correlation function, it may have lat- 
tice artifacts of order a and cannot be used to 
formulate an improvement condition. 

In terms of the bare unimproved correlation 
functions, 



(5) 
(6) 



as well as the bare improved correlation function. 



k\p^{xo,yQ) 
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e'^'^^{{A,)l{x){A^)l{y)Ql)X7) 



eq. may be rewritten, for ti < xq < as 

Zv[A:v(a;o) 4- acv9ofcT(a;o)] = 

Zl[k\K{t2. xo) - k\j^{tuXo)] + 0(a2) . (8) 

At zero quark mass m, this equation can be solved 
for cv- However, due to a peculiarity of the 
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Figure 1. Dependence of Yy on the quark mass. 



quenched approximation, i.e. the appearance of 
zero-modes in the quark propagator, it is not al- 
ways possible to simulate directly at zero quark 
mass For the particular boundary conditions 
and lattice size employed, here, this phenomenon 
is relevant when j3 — 6/.gg < 6.4. In this re- 
gion we have to perform an extrapolation to zero 
quark mass. However, in the course of the numer- 
ical simulations it turned out that cy, implicitly 
defined by eq. (H) also for finite to, is a steep func- 
tion of TO and could not be reliably extrapolated 
to zero quark mass when we were not able to sim- 
ulate at very small masses. The reason for this 
strong dependence is related to the fact that in 
the Ward identity (^) the mass term was left out: 
as a consequence, eq. is not even valid in the 
continuum limit. 

It is therefore natural to repeat the calculation, 
keeping the mass term. Eq. is then modified 
by an additional term and has corrections of order 
0(am) instead of 0(to). We solve the resulting 
equation formally for cy and denote the solution 
at finite mass by ly, which coincides with cy as 
TO ^ 0. As expected, this limit is now reached 
very smoothly (see figure |l|) . 

In a non-perturbative calculation, cy depends 
on the choices made for the various kinemat- 
ical parameters such as T/L or the particular 
Schrodinger Functional background field. Differ- 
ent choices lead to a variation of cy itself that is 
0(a). This is unavoidable and one must make a 



3 



specific choice; other choices would merely change 
the size of the O(a^) effects after improvement. 
On the other hand, one should search for val- 
ues of the kinematical parameters such that the 
O(a^) terms are small compared to the improve- 
ment term adokT^xo). To achieve this, we first 
of all choose ti and t2 in such a way that the 
time separations between the various fields in the 
correlation functions are as large as possible: 



ti = T/4, t2 = 3T/4 . 



(9) 



To fix the other parameters, we studied and tried 
to minimize the cutoff effects in eq. (^) at tree 
level of perturbation theory (the interested reader 
will find details in ref. pT[). 

With our final choices, we observe reasonably 
weak cutoff effects at tree level. Apart from this 
perturbative study, we also performed a rough nu- 
merical investigation at the largest values of the 
gauge coupling considering different improvement 
conditions (i.e. different choices of L/a and other 
parameters), and observed the extracted values 
for cv to be the same within errors. The final 
definition of cy includes an average over the three 
central time slices, xq, which reduces statistical 
uncertainties. 

3. RESULTS 

Our preliminary results are summarized in fig- 
ure 1^. Here we show the non-perturbatively com- 
puted values of cy as a function of the bare cou- 
pling (7q. The errors on cy receive a significant 
contribution from the (statistically independent) 
error on the ratio Zp^/Zy. 

In the region /3 < 6.4 we observe sizeable val- 
ues for Cy, far from 1-loop perturbation theory. 
In particular, at go = 1 the effect of the improve- 
ment term on the value of the vector meson de- 
cay constant can be as large as 10%. When 
approaches 0, we obtain agreement with the per- 
turbative result. 

This work is part of the ALPHA collaboration 
research programme. We thank DESY for allo- 
cating computer time for this project. 
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Figure 2. cy versus f/g. The dotted line is the 
l~loop perturbative result. 
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